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We  report  preparation  in  the  ground  state  of  collective  modes  of  motion  of  two  trapped  9Be+ 
ions.  This  is  a  crucial  step  towards  realizing  quantum  logic  gates  which  can  entangle  the  ions’ 
internal  electronic  states.  We  find  that  heating  of  the  modes  of  relative  ion  motion  is  substantially 
suppressed  relative  to  that  of  the  center-of-mass  modes,  suggesting  the  importance  of  these  modes  in 
future  experiments.  [S003 1-9007(98)06838-0] 


PACS  numbers:  03.67.Lx,  03.65. -w,  32.80.Pj 

In  physics,  quantum  computation  [1]  provides  a  gen¬ 
eral  framework  for  fundamental  investigations  into  sub¬ 
jects  such  as  entanglement,  quantum  measurement,  and 
quantum  information  theory.  Since  quantum  computation 
relies  on  entanglement  between  qubits,  any  implementa¬ 
tion  of  a  quantum  computer  must  offer  isolation  from 
the  effects  of  decoherence,  but  also  allow  controllable 
and  coherent  interaction  between  the  qubits.  Cirac  and 
Zoller  [2]  have  proposed  an  attractive  scheme  for  realiz¬ 
ing  a  quantum  computer,  which  is  scalable  to  an  arbitrary 
number  of  qubits.  Their  scheme  is  based  on  a  collection 
of  trapped  atomic  ions,  where  each  qubit  (one  per  ion) 
is  comprised  of  a  pair  of  the  ions’  internal  states,  while 
quantum  information  is  transferred  between  different  ions 
using  a  particular  quantized  mode  of  the  ions’  collective 
motion.  This  “quantum  data  bus”  must  first  be  initialized 
in  a  pure  quantum  state  [2]:  for  example,  its  ground  state 
[3].  The  basics  of  this  scheme  have  been  demonstrated 
experimentally  in  a  fundamental  logic  gate  (a  Controlled- 
NOT)  operating  between  a  motional  mode  of  a  single 
trapped  ion  and  two  of  the  ion’s  internal  states  [4].  In 
that  work,  the  motional  state  was  initialized  in  the  ground 
state  by  laser  cooling  [5].  The  next  step  towards  imple¬ 
menting  the  Cirac-Zoller  scheme  is  to  cool  at  least  one 
mode  of  collective  motion  of  multiple  ions  to  the  ground 
state.  In  this  Letter,  we  describe  the  first  experiments  to 
realize  this  goal.  We  also  report  a  significant  difference 
between  the  decoherence  rates  of  the  center-of-mass  and 
non-center-of-mass  modes  of  motion. 

We  confine  9Be+  ions  in  a  coaxial-resonator-based  rf 
(Paul)  trap,  similar  to  that  described  in  Ref.  [6].  The 
electrodes  in  this  trap  are  made  from  125 -yam- thick 
sheets  of  Be  metal,  as  shown  in  Fig.  1.  We  apply  a  po¬ 
tential  cf)(t )  =  Vocos(ft7T)  +  U0  to  the  (elliptical)  ring 
electrode  relative  to  the  end  cap  electrodes.  If  several 
ions  are  trapped  and  cooled,  they  will  naturally  align 
themselves  along  the  major  axis  of  the  ring  electrode. 
The  electrode’s  elliptical  shape,  in  combination  with 
U{)  >  0,  allows  a  linear  crystal  to  be  maintained  while 
suppressing  rf  micromotion  of  the  ions  along  this  direc¬ 
tion  [7].  With  V0  «  520  V,  CIt/2tt  -  238  MHz,  and 
Uq  =  0  V,  the  pseudopotential  oscillation  frequencies  are 


(0)x,  coy,  coz)/2tt  «  (4.6, 12.7, 17.0)  MHz.  With  U0  = 
18.2  V,  the  frequencies  become  (8.6, 17.6, 9.3)  MHz. 
Figure  1  shows  two  ions  confined  in  the  trap  and  imaged 
with  an  f/3  lens  system  onto  a  position-sensitive  photo¬ 
multiplier  tube. 

The  ions  are  cooled  and  probed  with  laser  beams  whose 
geometry  is  indicated  in  Fig.  2(a).  The  relevant  level 
structure  of  9Be+  is  shown  in  Fig.  2(b).  The  quan¬ 
tization  axis  is  defined  by  an  applied  static  magnetic 
field;  \B\  —  0.2  mT.  The  levels  of  interest  for  quan¬ 
tum  logic  operations  are  the  2s  2S\/2\F  =  2,  m?  =  2)  and 
2s2S\/2\F  =  1  =  1)  states,  abbreviated  by  1 1)  and 

It),  respectively.  Laser  beams  Dl,  D2,  and  D3  are  a+ 
polarized  and  focused  to  nearly  saturate  the  ions  (7sat  ~ 
85  mW cm-2).  Beams  Dl  and  D2  provide  Doppler  pre- 
cooling  in  all  three  dimensions,  and  beam  D3  prevents  op¬ 
tical  pumping  to  the  |F  =  2,  m?  =  1)  state.  The  ||)  — ► 
2p  ^7*3/2 1 E  =  3 =  3)  transition  (radiative  linewidth 
y I 2tt  —  19.4  MHz),  driven  by  D2,  is  a  cycling  transi¬ 
tion,  which  allows  us  to  detect  the  ion’s  electronic  state 
(| |)  or  It))  with  nearly  unit  detection  efficiency. 

Beams  R1  (a+/a~  polarized)  and  R2  (7 t  polarized) 
are  used  to  drive  stimulated  Raman  transitions  between 


1 — i 

3  pm 


FIG.  1  (color).  Two  ions  trapped  in  an  elliptical  rf  (Paul)  trap. 
The  ring  has  an  aspect  ratio  of  3:2  and  the  major  axis  is 
525  /am  long.  The  slot  which  forms  the  end  caps  is  250  /am 
across.  A  potential  <p(t)  is  applied  to  the  ring  (see  text).  The 
Be  sheets  are  —125  /am  thick.  With  an  x-axis  pseudopotential 
oscillation  frequency  coxI2tt  =  4.6  MHz,  the  ion-ion  spacing 
is  approximately  3  /am. 
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1 1)  and  1 1),  through  the  virtual  2 p  2P\/2  state  [5].  These 
beams  are  derived  from  a  single  laser,  whose  output  is 
split  by  an  acousto-optic  modulator  [8].  The  beams  are 
detuned  by  \/2ir  —  40  GHz  to  the  red  of  the  2s  2 S 1/2  — ► 
2p2P\/2  transition,  and  their  frequency  difference  is 
tuned  around  the  2s  2S 1/2  hyperfine  splitting  of  coo/2tt  — 
1.25  GHz.  (Here,  c 00  includes  stable  shifts  of  a  few  mega¬ 
hertz  from  the  Zeeman  and  ac  Stark  effects.)  R2  is  di¬ 
rected  along  (-1/V2)*  +  (1/2)  (— y  +  z).  If  R1  _L  R2 
as  in  Fig.  2,  then  the  Raman  beam  wave  vector  difference 
8k  ||  x ,  and  the  transitions  are  sensitive  to  ion  motion 
only  in  this  direction.  If,  however,  R1  is  counterpropa- 
gating  to  R2,  the  transitions  become  sensitive  to  motion 
in  all  three  dimensions. 

When  two  cold  ions  are  held  in  the  trap  and  undergo 
small  oscillations  about  their  equilibrium  positions,  we 
may  solve  the  equations  of  motion  using  normal  mode 
coordinates.  For  two  ions  lying  along  the  x  axis  there  are 
two  modes  involving  motion  along  this  axis:  the  center-of- 
mass  (COM)  mode  (in  which  the  ions  move  together  with 
frequency  &>com  =  cox)  and  the  stretch  mode  (wherein 
the  ions  move  out  of  phase,  with  frequency  costY  = 
V^^com)*  The  other  motional  frequencies  are  coy  (y 
center  of  mass),  coz  (z  center  of  mass),  j(o2  —  co2  (xy 
rocking),  and  s/(o2z  —  a>2  (xz  rocking). 

The  lower  traces  in  Fig.  3,  taken  with  8k  ||  x,  show 
an  x-axis  normal  mode  spectrum;  results  for  the  y  and 
z  modes  are  very  similar.  We  take  the  data  with  the 
following  steps:  first  we  turn  on  beams  Dl,  D2,  and 
D3  for  approximately  10  yccs  to  Doppler  cool  the  ions 
to  the  Lamb-Dicke  regime,  where  the  ions’  confinement 
is  much  smaller  than  the  laser  wavelength.  Next,  we 


turn  off  beam  D2,  and  leave  beams  Dl  and  D3  on  for 
5  juls  to  optically  pump  both  ions  to  the  |  j)  state.  We 
then  turn  on  only  the  Raman  beams  R1  and  R2  for  a 
time  fpr,  with  relative  detuning  coq  +  <5pr  (the  “Raman 
probe”  pulse).  Finally,  we  drive  the  cycling  transition 
with  D2  and  measure  the  ions’  fluorescence.  We  repeat 
the  experiment  at  a  rate  of  a  few  kilohertz  while  slowly 
sweeping  <5pr.  If  the  Raman  beam  difference  frequency 
is  resonant  with  a  transition,  then  an  ion  is  driven  from 
|  j)  — ►  It)  and  the  D2-driven  fluorescence  rate  drops. 

For  a  single  ion,  the  carrier  transition  (<5pr  =  0)  causes 
the  population  to  undergo  sinusoidal  Rabi  oscillations 
between  |j)  and  ||)  [9].  The  effective  Rabi  frequency 
is  12  =  gig2/A  —  27 t  X  250  kHz,  where  g  1,  g2  are  the 
single-photon  resonant  Rabi  frequencies  of  beams  R1  and 
R2.  (We  assume  A  »  y,  com  »  H,  where  com  is  the  fre¬ 
quency  of  the  motional  mode  of  interest.)  If  <5pr  =  —  cox 
(the  first  lower  x  sideband),  then  the  transition  couples  the 
states  || ,nx)  and  \\,nx  —  1),  where  nx  is  the  vibrational 
level  of  the  quantized  motion  along  x.  In  the  Lamb- 
Dicke  regime,  the  corresponding  Rabi  frequency  is  given 
by  £lnx,nx- 1  =  VxyJrhQ  [9].  Here,  r)x  =  x0\8k  •  x\ 
is  the  Lamb-Dicke  parameter  (=  0.23  when  o)xI2tt  = 
8.6  MHz)  and  xo  =  V H/(2mcox )  is  the  spread  of  the 
nx  =  0  wave  function,  with  m  being  the  ion’s  mass). 
(Note  that  if  the  ion  is  in  the  nx  =  0  state  of  motion,  this 
lower  sideband  vanishes.)  The  first  upper  x  sideband 
transition  (<5pr  =  +cox)  couples  1 1,  nx)  and  ||,  nx  +  1) 
with  Rabi  frequency  Plnx,nx+\  =  Vxs/nx  +  1  fT 

In  the  case  of  two  ions  driven  on  the  carrier  transi¬ 
tion,  each  ion  independently  undergoes  Rabi  oscillations 


FIG.  2.  (a)  Laser  beam  geometry.  The  trap  ring  electrode  is 

shown  rotated  45°  out  of  the  page.  The  end  cap  electrodes 
are  omitted  for  clarity  (see  Fig.  1).  A  magnetic  field  B 
of  magnitude  0.2  mT  defines  the  quantization  axis  along 
—  (1/V2)x  +  (1/2)  (y  —  z),  and  laser  beam  polarizations  are 
indicated,  (b)  Relevant  9Be+  energy  levels  (not  to  scale), 
indicated  by  F,  mF  quantum  numbers  in  the  ground  state.  2P 
fine  structure  splitting  is  —197  GHz,  2  S 1/2  hyperfine  splitting  is 
coo/2tt  —  1.25  GHz,  2Pi/2  hyperfine  splitting  is  —237  MHz, 
and  the  2P^/2  hyperfine  structure  («y/27r  —  19.4  MHz)  is 
not  resolved.  All  optical  transitions  are  near  A  —  313  nm,  and 
A/2t t  «  40  GHz. 
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FIG.  3.  Spectrum  of  sidebands  due  to  two-ion  x-axis  normal 
mode  motion:  (from  left  to  right)  lower  stretch,  lower  COM, 
upper  COM,  and  upper  stretch.  The  ordinate  is  the  detuning 
of  the  Raman  probe  beam  difference  frequency  from  the 
carrier  transition.  The  abscissa  shows  the  ion  fluorescence 
(proportional  to  the  expectation  value  of  the  number  of  atoms  in 
the  state  1 1)),  plus  a  constant  background  (whose  approximate 
level  for  the  lower  curves  is  indicated  by  the  dashed  line). 
The  solid  lines,  meant  as  guides  to  the  eye,  are  fits  to 
Gaussians.  The  lower  traces  show  the  effects  of  Doppler 
cooling.  The  upper  traces,  offset  vertically  for  clarity,  show  the 
effects  of  several  pulses  of  Raman  cooling  on  the  mode  which 
is  displayed.  Vanishing  lower  motional  sidebands  indicate 
cooling  to  the  ground  state  of  motion.  The  peak  widths  are 
consistent  with  the  Raman  probe  pulse  lengths  (—3  /ms). 
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between  |  j)  and  |  |)  with  Rabi  frequency  II .  Since 
the  laser  beam  waists  (—20  ^m)  are  much  larger 
than  the  ion-ion  separation  (—2  jmm),  the  ions  are 
equally  illuminated.  Nonetheless,  if  the  micromotion 
of  the  two  ions  is  different,  then  the  reduction  of  the 
carrier  (and  sideband)  transition  strengths  due  to  the 
micromotion  will  give  a  different  Rabi  frequency  for 
each  ion  [7,10].  This  could  be  used  as  a  means  of 
selectively  addressing  the  ions  [11];  however,  in  the 
present  work  the  two  ions’  Rabi  frequencies  were 
equal. 


Since  the  sideband  transitions  affect  the  motional  state, 
which  is  a  shared  property  of  both  ions,  such  transitions 
produce  entanglement  between  the  ions’  spins  and  their 
collective  motion  [12].  The  system  can  no  longer  be 
treated  as  two,  independent,  two-level  systems  and  the 
measured  fluorescence  following  the  Raman  probe  is  a 
complicated  function  of  the  probe  pulse  duration  £pr.  For 
example,  given  an  initial  state  (where  n  is  the 

vibrational  level  of  the  COM  or  stretch  motion  along  the 
v  axis)  driven  on  the  corresponding  lower  sideband  for  a 
time  fpr,  the  wave  function  evolves  as 


\^n{tm))  =  1 


2  n  -  1 


[1  -  COs(G?nr)] 


II,  In) 


ie 


2  n 
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sin(Gfpr) 


(IU)  ±  <IU))k  ~  1) 
V2 


+  e 


ie 


■Jn2  —  n 
2n  —  1 


[1 


cos(Gfpr)]|T,t,«  “  2), 
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where  G  =  ^2(2 n  —  1)  fl  r\x,m  and  0,  (f)  are  the  sum  and 
difference  of  the  Raman  beam  phases  at  the  ions.  On  the 
COM  sideband  [top  sign  in  Eq.  (1)],  r\x  m  =  77x,com  = 
r]x/yj2  (down  by  a  factor  of  \j2  from  the  single-ion  case 
due  to  the  extra  mass  of  the  two-ion  string),  whereas  on  the 
stretch  sideband  (lower  sign),  r]x  m  =  r]x  stY  =  rjx/s! 2^/3. 
The  expressions  for  transitions  on  the  upper  motional 
sidebands  are  similar.  If,  before  the  Raman  probe  pulse, 
the  ions  have  probability  pn  of  being  in  the  motional  state 
| ft),  the  subsequently  measured  average  fluorescence  from 
the  cycling  transition  is 

S(tpr)  =  X  Pn(2\(ll,n\tffn(tpl))\2 
n 

+  Kl.t.n  -  l|*An0pr)>|2 

+  KU,n  -  1|«A»M>I2)-  (2) 

This  signal  is  proportional  to  the  expectation  value  of  the 
number  of  atoms  in  the  state  |  j).  For  the  data  shown  in 
Fig.  3,  tVY  was  chosen  to  maximize  the  sideband  features. 

The  upper  traces  in  Fig.  3  show  the  effects  of  adding 
several  cycles  of  Raman  cooling  [5]  on  one  particular 
v  mode  after  the  Doppler  cooling  but  before  the  probe 
pulse.  The  reduction  in  the  mean  vibrational  number  (n) 
is  indicated  by  the  reduction  in  size  of  the  lower  sideband, 
which  vanishes  in  the  limit  ( n )  — ►  0.  The  data  are 
consistent  with  a  thermal  state  of  (^com)  =  0*H-ao3  or 
(nstr)  =  O.Olloiou  This  implies  that  the  COM  and  stretch 

T  3  T 1 

modes  are  in  their  ground  states  90  -12%  and  99-7%  of  the 
time,  respectively.  We  have  also  simultaneously  cooled 
the  COM  and  stretch  modes  along  x,  to  comparable  values 
of  (n)  (and  have  separately  cooled  the  other  four  motional 
modes — y  and  z  COM,  xy  rocking,  and  xz  rocking — to 
near  their  ground  states). 

Each  cycle  of  Raman  cooling  consists  of  (i)  a  pulse  of 
the  Raman  beams  with  their  difference  frequency  tuned  to 


one  of  the  lower  sidebands  (COM  or  stretch  mode)  and 
(ii)  optical  repumping  to  the  |  j)  state  driven  by  beams  D1 
and  D3.  The  Raman  transition  reduces  the  vibrational  en¬ 
ergy  by  Hcom,  whereas  the  repumping,  on  average,  heats 
each  mode  by  approximately  the  recoil  energy 
Therefore,  the  ion  is  cooled  through  this  process.  Five 
pulses  of  Raman  cooling  were  used  for  the  data  shown 
in  Fig.  3.  The  exact  durations  of  the  Raman  pulses  were 
chosen  to  optimize  the  cooling  rate — each  pulse  was  ap¬ 
proximately  5  juls  long. 

For  an  ion-trap  implementation  of  a  quantum  computer, 
the  motional  modes  are  most  susceptible  to  decoherence. 
The  ions’  motional  states  lose  coherence  if  they  couple  to 
(stochastic)  electric  fields  caused  by  fluctuating  potentials 
on  the  electrodes.  This  leads  to  heating,  which  has 
previously  been  observed  in  single  ions  [5,10,13];  in 
Ref.  [5],  the  heating  drove  the  ion  out  of  the  motional 
(COM)  ground  state  in  approximately  1  ms.  We  have 
performed  similar  heating  measurements  on  the  COM  and 
non-COM  modes  of  motion  of  two  ions.  The  results  are 
summarized  in  Table  I.  The  heating  rate  was  determined 
by  inserting  a  delay  between  laser  cooling  and  the  Raman 

TABLE  I.  Heating  rates  of  the  six  normal  modes  of  two 
trapped  ions.  The  Raman  beams  were  counterpropagating  for 
the  y-  and  z-axis  data,  making  the  Raman  probe  sensitive  to 
motion  in  all  three  dimensions.  Note  that  the  COM  modes  are 
heated  at  a  much  higher  rate  than  the  non-COM  modes  (see 
text).  (The  precision  with  which  the  heating  rates  are  given  for 
the  last  five  modes  is  limited  by  measurement  noise.) 


Mode 

&>m/27r  (MHz) 

8(n)/8t  (ms  !) 

vcom 

8.6 

19^13 

ycoM 

17.6 

>10 

ZCOM 

9.3 

>20 

Vstr 

14.9 

<0.18 

•Vyrocking 

15.4 

<1 

VZrocking 

3.6 

<0.5 
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probe.  The  main  results  from  these  data  are  that  the  COM 
modes  are  heated  out  of  the  ground  state  much  more 
quickly  than  the  non-COM  modes.  This  can  be  explained 
as  follows. 

The  COM  modes,  in  which  both  ions  move  in  phase, 
can  be  excited  by  a  uniform  electric  field.  However, 
no  non-COM  mode  can  be  excited  by  a  uniform  electric 
field  [14] — since  these  modes  involve  differential  motion 
of  the  ions,  they  can  be  driven  only  by  field  gradients. 
If  the  fluctuating  field  at  the  ion  (along  the  direction 
of  motion  of  the  mode  of  interest)  is  E(t ),  an  estimate 
of  the  corresponding  field  gradient  is  E{t)/d,  where  d 
is  a  characteristic  internal  dimension  of  the  trap.  For 
stochastic  fields,  the  COM  heating  rate  scales  as  (E2(t))\ 
the  non-COM  mode  heating  rates  scale  as  ([^pAv]2) 
(where  Av  is  the  ion-ion  separation),  down  by  a  factor 
of  104  for  the  present  trap.  Similarly,  other  non-COM 
modes  for  more  than  two  ions  can  be  excited  only  by 
higher-order  field  gradients,  leading  to  further  reductions 
in  their  heating. 

This  suggests  using  non-COM  modes  for  the  quantum 
data  bus  in  the  Cirac-Zoller  scheme.  Excitation  of  the 
“spectator”  COM  modes  along  the  direction  to  which  the 
Raman  transitions  are  sensitive  will  still  alter  the  Rabi 
frequencies,  but  these  effects  will  be  higher  order  in  the 
Lamb-Dicke  parameter  [10].  In  the  two-ion  example,  in 
the  Lamb-Dicke  regime,  the  Rabi  frequency  for  a  first 
sideband  transition  \n\)  — ►  \n[ )  on  (cold)  mode  1,  given 
that  (hotter)  mode  2  is  in  the  state  Inf),  is  [10] 

Qnun'M 2)  =  e~{v'+Vl)/2(l  -  n 2yl) ,  (3) 

where  n\>  denotes  the  larger  of  n\  or  n\,  and  i)\  and 
772  are  the  Lamb-Dicke  parameters  for  modes  1  and  2, 
respectively.  Fluctuations  in  the  Rabi  frequency  of  mode 
1  due  to  fluctuations  in  n2  therefore  occur  in  order  r\\. 
However,  for  the  conditions  of  the  present  experiment, 
even  if  quantum  logic  operations  were  performed  using 
the  v-stretch  mode,  the  x-COM  mode  heating  would  still 
limit  the  number  of  operations  to  around  ten  by  the  above 
mechanism.  Clearly,  this  heating  must  be  eliminated  in 
future  experiments. 

The  two-ion  cooling  results  presented  here  are  compa¬ 
rable  to  our  previous  single-ion  results  [5],  indicating  that 
rf  heating  should  not  be  a  concern  for  small  numbers  of 
ions  [10].  Comparable  cooling  for  N  >  2  ions  should  not 
present  any  fundamental  difficulties,  as  long  as  spurious 
overlaps  of  motional  modes  are  avoided. 

The  preparation  of  a  pure  state  of  motion  (the  ground 
state)  of  multiple  trapped  ions  represents  the  first  step 
towards  realizing  quantum  logic  operations  on  them. 


Such  operations  should  lead  to  the  creation  of  arbitrary 
entangled  states  of  massive  particles,  including  EPR-  or 
GHZ-like  spin  states  [15].  Unlike  other  systems  in  which 
EPR  states  have  been  generated,  it  should  be  possible  to 
reliably  create  these  states  on  demand  [11]  rather  than  by 
a  selection  process,  and  to  detect  them  with  nearly  perfect 
efficiency  [16]. 
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